The impact of small movements with dual lumen cannulae during venovenous extracorporeal membrane oxygenation: A computational fluid dynamics analysis.

BACKGROUND AND OBJECTIVES: Venovenous Extracorporeal Membrane Oxygenation (VV ECMO) provides respiratory support to patients with severe lung disease failing conventional medical therapy. An essential component of the ECMO circuit are the cannulas, which drain and return blood into the body. Despite being anchored to the patient to prevent accidental removal, minor cannula movements are common during ECMO. The clinical and haemodynamic consequences of these small movements are currently unclear. This study investigated the risk of thrombosis and recirculation caused by small movements of a dual lumen cannula (DLC) in an adult using computational fluid dynamics. METHODS: The 3D model of an AVALON Elite DLC (27 Fr) and a patient-specific vena cava and right atrium were generated for an adult patient on ECMO. The baseline cannula position was generated where the return jet enters the tricuspid valve. Alternative cannula positions were obtained by shifting the cannula 5 and 15 mm towards inferior (IVC) and superior (SVC) vena cava, respectively. ECMO settings of 4 L/min blood flow and pulsatile flow at SVC and IVC were applied. Recirculation was defined as a scalar value indicating the infused oxygenated blood inside the drainage lumen, while thrombosis risk was evaluated by shear stress, stagnation volume, washout, and turbulent kinetic energy. RESULTS: Recirculation for all models was less than 3.1 %. DLC movements between -5 to 15 mm increased shear stress and turbulence kinetic energy up to 24.7 % and 11.8 %, respectively, compared to the baseline cannula position leading to a higher predicted thrombosis risk. All models obtained a complete washout after nine seconds except for when the cannula migrated 15 mm into the SVC, indicating persisting stasis and circulating zones. CONCLUSION: In conclusion, small DLC movements were not associated with an increased risk of recirculation. However, they may increase the risk of thrombosis due to increased shear rate, turbulence, and slower washout of blood. Developing effective cannula securement devices may reduce this risk.

Blood flow and emboli transport patterns during venoarterial extracorporeal membrane oxygenation: A computational fluid dynamics study.

PROBLEM: Despite advances in Venoarterial Extracorporeal Membrane Oxygenation (VA-ECMO), a significant mortality rate persists due to complications. The non-physiological blood flow dynamics of VA-ECMO may lead to neurological complications and organ ischemia. Continuous retrograde high-flow oxygenated blood enters through a return cannula placed in the femoral artery which opposes the pulsatile deoxygenated blood ejected by the left ventricle (LV), which impacts upper body oxygenation and subsequent hyperoxemia. The complications underscore the critical need to comprehend the impact of VA-ECMO support level and return cannula size, as mortality remains a significant concern. AIM: The aim of this study is to predict and provide insights into the complications associated with VA-ECMO using computational fluid dynamics (CFD) simulations. These complications will be assessed by characterising blood flow and emboli transport patterns through a comprehensive analysis of the influence of VA-ECMO support levels and arterial return cannula sizes. METHODS: Patient-specific 3D aortic and major branch models, derived from a male patient's CT scan during VA-ECMO undergoing respiratory dysfunction, were analyzed using CFD. The investigation employed species transport and discrete particle tracking models to study ECMO blood (oxygenated) mixing with LV blood (deoxygenated) and to trace emboli transport patterns from potential sources (circuit, LV, and aorta wall). Two cannula sizes (15 Fr and 19 Fr) were tested alongside varying ECMO pump flow rates (50%, 70%, and 90% of the total cardiac output). RESULTS: Cannula size did not significantly affect oxygen transport. At 90% VA-ECMO support, all arteries distal to the aortic arch achieved 100% oxygen saturation. As support level decreased, oxygen transport to the upper body also decreased to a minimum saturation of 73%. Emboli transport varied substantially between emboli origin and VAECMO support level, with the highest risk of cerebral emboli coming from the LV with a 15 Fr cannula at 90% support. CONCLUSION: Arterial return cannula sizing minimally impacted blood oxygen distribution; however, it did influence the distribution of emboli released from the circuit and aortic wall. Notably, it was the support level alone that significantly affected the mixing zone of VA-ECMO and cardiac blood, subsequently influencing the risk of embolization of the cardiogenic source and oxygenation levels across various arterial branches.

Effect of RVAD Cannulation Length on Right Ventricular Thrombosis Risk: An In Silico Investigation.

Left ventricular assist devices (LVADs) have been used off-label as long-term support of the right heart due to the lack of a clinically approved durable right VAD (RVAD). Whilst various techniques to reduce RVAD inflow cannula protrusion have been described, the implication of the protrusion length on right heart blood flow and subsequent risk of thrombosis remains poorly understood. This study investigates the influence of RVAD diaphragmatic cannulation length on right ventricular thrombosis risk using a patient-specific right ventricle in silico model validated with particle image velocimetry. Four cannulation lengths (5, 10, 15 and 25 mm) were evaluated in a one-way fluid-structure interaction simulation with boundary conditions generated from a lumped parameter model, simulating a biventricular supported condition. Simulation results demonstrated that the 25-mm cannulation length exhibited a lower thrombosis risk compared to 5-, 10- and 15-mm cannulation lengths due to improved flow energy distribution (25.2%, 24.4% and 17.8% increased), reduced stagnation volume (72%, 68% and 49% reduction), better washout rate (13.0%, 11.6% and 9.1% faster) and lower blood residence time (6% reduction). In the simulated scenario, our findings suggest that a longer RVAD diaphragmatic cannulation length may be beneficial in lowering thrombosis risk; however, further clinical studies are warranted.

Flow characterization of Maquet and Bio-Medicus multi-stage drainage cannulae during venoarterial extracorporeal membrane oxygenation.

BACKGROUND: Drainage cannulae extract blood from a patient during venoarterial extracorporeal membrane oxygenation (VA ECMO), a treatment that temporarily supports patients undergoing severe heart and/or lung dysfunction. Currently, the two most commonly used multi-stage drainage cannulae are manufactured by Maquet and Bio-Medicus, but their designs vary in many aspects which impacts the generated flow dynamics. Therefore, this study aimed to use computational fluid dynamics (CFD) to explore the flow characteristics of the aforementioned cannulae and their impact on complications such as thrombosis. METHODS: The Maquet and Bio-Medicus cannulae were 3D modelled within a patient-specific geometry of the venous vasculature taken from a computed tomography scan of a patient undergoing VA ECMO. A drainage flow rate of 4 L/min was assigned to each cannula. Lastly, a stress blended eddy simulation turbulence model was employed to resolve bulk flow turbulence. RESULTS: The proximal row of side holes in both cannulae generated high intensity counter-rotating vortices, thus generating supraphysiological shear. These proximal rows were also responsible for the majority of flow extraction in both cannulae (>1.6 L/min). Despite identical simulation settings, each cannulae had differing impacts on global flow dynamics. For instance, the Bio-Medicus model produced a total stagnant blood volume of 25.6 ml, compared to 17.8 ml the Maquet cannula, thereby increasing the risk of thrombosis. CONCLUSIONS: Overall, our results demonstrate that differences in design clearly impact flow dynamics and risk of complications. Therefore, further work in optimizing cannula design may be beneficial to prevent harmful flow characteristics.

Computational analysis of thrombosis risk with variations in left ventricular assist device inflow cannula design in a multi-patient model.

BACKGROUND AND OBJECTIVES: Left ventricular assist devices (LVADs) are mechanical pumps used to support patients with end-stage heart failure. The inflow cannula is a critical component of the LVAD as it connects the pump to the left ventricle, allowing blood to be drawn from the heart. However, the design of the cannula can significantly impact LV hemodynamics and cause complications, including thrombosis. Therefore, this study aimed to analyze the numerical effects of left ventricle (LV) size on cannula design in order to enhance hemodynamic performance using post-operative left ventricular assist device (LVAD) models. METHODS: A parametric design evaluation of two different inflow cannulas were carried out on left ventricles (LV) of varying sizes (ranging from 154 to 430 ml) constructed from computerized tomography (CT) data from VAD patients using computational fluid dynamics (CFD) simulations. The study analyzed three key factors contributing to thrombosis formation: blood residence time, blood stagnation ratio, and wall shear stress. RESULTS: Results showed higher blood residence time and stagnation ratio for larger left ventricular sizes. In addition, increasing the insertion length of the cannula reduced the average wall shear stress. CONCLUSION: Overall, the study's findings suggest that the optimal cannula shape for LVADs varies with left ventricular size.

The effect of drainage cannula tip position on risk of thrombosis during venoarterial extracorporeal membrane oxygenation.

BACKGROUND AND OBJECTIVES: Venoarterial extracorporeal membrane oxygenation (VA ECMO) is able to support critically ill patients undergoing refractory cardiopulmonary failure. It relies on drainage cannulae to extract venous blood from the patient, but cannula features and tip position may impact flow dynamics and thrombosis risk. Therefore, this study aimed to investigate the effect of tip position of single-stage (SS) and multi-stage (MS) VA ECMO drainage cannulae on the risk of thrombosis. METHODS: Computational fluid dynamics was used to model flow dynamics within patient-specific geometry of the venous vasculature. The tip of the SS and MS cannula was placed in the superior vena cava (SVC), SVC-Right atrium (RA) junction, mid-RA, inferior vena cava (IVC)-RA junction, and IVC. The risk of thrombosis was assessed by measuring several factors. Blood residence time was measured via an Eulerian approach through the use of a scalar source term. Regions of stagnant volume were recognised by identifying regions of low fluid velocity and shear rate. Rate of blood washout was calculated by patching the domain with a scalar value and measuring the rate of fluid displacement. Lastly, wall shear stress values were determined to provide a qualitative understanding of potential blood trauma. RESULTS: Thrombosis risk varied substantially with position changes of the SS cannula, which was less evident with the MS cannula. The SS cannula showed reduced thrombosis risk arising from stagnant regions when placed in the SVC or SVC-RA junction, whereas an MS cannula was predicted to create stagnant regions during all tip positions. When positioned in the IVC-RA junction or IVC, the risk of thrombosis was higher in the SS cannula than in the MS cannula due to both high and low shear flow. CONCLUSION: Tip position of the drainage cannula impacts cannula flow dynamics and, subsequently, the risk of thrombosis. The use of MS cannulae can reduce high shear-related thrombosis, but SS cannulae can eliminate stagnant regions when advanced into the SVC. Therefore, the choice of cannula design and tip position should be carefully considered during cannulation.

The effect of arterial cannula tip position on differential hypoxemia during venoarterial extracorporeal membrane oxygenation.

Interaction between native ventricular output and venoarterial extracorporeal membrane oxygenation (VA ECMO) flow may hinder oxygenated blood flow to the aortic arch branches, resulting in differential hypoxemia. Typically, the arterial cannula tip is placed in the iliac artery or abdominal aorta. However, the hemodynamics of a more proximal arterial cannula tip have not been studied before. This study investigated the effect of arterial cannula tip position on VA ECMO blood flow to the upper extremities using computational fluid dynamics simulations. Four arterial cannula tip positions (P1. common iliac, P2. abdominal aorta, P3. descending aorta and P4. aortic arch) were compared with different degrees of cardiac dysfunction and VA ECMO support (50%, 80% and 90% support). P4 was able to supply oxygenated blood to the arch vessels at all support levels, while P1 to P3 only supplied the arch vessels during the highest level (90%) of VA ECMO support. Even during the highest level of support, P1 to P3 could only provide oxygenated VA-ECMO flow at 0.11 L/min to the brachiocephalic artery, compared with 0.5 L/min at P4. This study suggests that cerebral perfusion of VA ECMO flow can be increased by advancing the arterial cannula tip towards the aortic arch.

Computational investigation of nasal surface coverage from squeeze bottle and Neti Pot saline irrigation flow.

BACKGROUND AND OBJECTIVE: Nasal saline irrigation is a common therapy for inflammatory nasal and paranasal disease or for managing post nasal and sinus surgery recovery. Two common irrigation devices include the netipot and squeeze bottles, where anecdotally, these devices alleviate congestion, facial pain, and pressure. However, a quantitative evaluation of these devices' performance and the fluid dynamics responsible for the irrigation distribution through the nose is lacking. This study tracked the liquid surface coverage and wall shear stresses during nasal saline irrigation produced from a Neti Pot and squeeze bottle. METHODS: This study used transient computational fluid dynamics (CFD) simulations to investigate the saline irrigation flow field in a subject-specific sinonasal model. The computational nasal cavity model was constructed from a high-resolution computed tomography scan (CT). The irrigation procedure applied a head position tilted at 90° forward using an 80 ml squeeze bottle and 120 ml Neti Pot. RESULTS: The results from a single sinonasal model demonstrated that the Neti Pot irrigation was more effective in delivering saline solution to the nasal cavity on the contralateral side of irrigation due to typically larger volumes but at the expense of reduced flow and shearing rates, as the flow entered under gravitational forces. The squeeze bottle irrigation provided greater surface coverage on the side of irrigation. CONCLUSIONS: The results from the single patient model, demonstrated the Neti Pot increased surface coverage in the paranasal sinuses. Reducing the jet diameter may aid the direct targeting of a specific region at the side of irrigation by preventing the impingement of the jet to the nasal passage surface and redirection of the flow. Evaluating this performance across a wider cohort of patients can strengthen the findings.

Nasal Irrigation Delivery in Three Post-FESS Models From a Squeeze-bottle Using CFD.

PURPOSE: Nasal saline irrigation is highly recommended in patients following functional endoscopic sinus surgery (FESS) to aid the postoperative recovery. Post-FESS patients have significantly altered anatomy leading to markedly different flow dynamics from those found in pre-op or non-diseased airways, resulting in unknown flow dynamics. METHODS: This work investigated how the liquid stream disperses through altered nasal cavities following surgery using Computational Fluid Dynamics (CFD). A realistic squeeze profile was determined from physical experiments with a 27-year-old male using a squeeze bottle with load sensors. The administration technique involved a head tilt of 45-degrees forward to represent a head position over a sink. After the irrigation event that lasted 4.5 s, the simulation continued for an additional 1.5 s, with the head orientation returning to an upright position. RESULTS: The results demonstrated that a large maxillary sinus ostium on the right side allows saline penetration into this sinus. The increased volume of saline entering the maxillary sinus limits the saline volume available to the rest of the sinonasal cavity and reduces the surface coverage of the other paranasal sinuses. The average wall shear stress was higher on the right side than on the other side for two patients. The results also revealed that head position alters the sinuses' saline residual, especially the frontal sinuses. CONCLUSION: While greater access to sinuses is achieved through FESS surgery, patients without a nasal septum limits posterior sinus penetration due to the liquid crossing over to the contralateral cavity and exiting the nasal cavity early.

Wet surface wall model for latent heat exchange during evaporation.

Air conditioning is a dual heat and mass transfer process, and the human nasal cavity achieves this through the mucosal wall surface, which is supplied with an energy source through the sub-epithelial network of capillaries. Computational studies of air conditioning in the nasal cavity have included temperature and humidity, but most studies solved these flow parameters separately, and in some cases, a constant mucosal surface temperature was used. Recent developments demonstrated that both heat and mass transfer need to be modeled. This work expands on existing modeling efforts in accounting for the nasal cavity's dual heat and mass transfer process by introducing a new subwall model, given in the Supplementary Materials. The model was applied to a pipe geometry, and a human nasal cavity was recreated from CT-scans, and six inhalation conditions were studied. The results showed that when the energy transfer from the latent heat of evaporation is included, there is a cooling effect on the mucosal surface temperature.

N95 respirator mask breathing leads to excessive carbon dioxide inhalation and reduced heat transfer in a human nasal cavity.

Face masks and respirators are used to filter inhaled air, which may contain airborne droplets and high particulate matter (PM) concentrations. The respirators act as a barrier to the inhaled and exhaled air, which may change the nasal airflow characteristics and air-conditioning function of the nose. This study aims to investigate the nasal airflow dynamics during respiration with and without an N95 respirator driven by airflow through the nasal cavity to assess the effect of the respirator on breathing conditions during respiration. To achieve the objective of this study, transient computational fluid dynamics simulations have been utilized. The nasal geometry was reconstructed from high-resolution Computed Tomography scans of a healthy 25-year-old female subject. The species transport method was used to analyze the airflow, temperature, carbon dioxide (CO2), moisture content (H2O), and temperature distribution within the nasal cavity with and without an N95 respirator during eight consecutive respiration cycles with a tidal volume of 500 ml. The results demonstrated that a respirator caused excessive CO2 inhalation by approximately 7 × greater per breath compared with normal breathing. Furthermore, heat and mass transfer in the nasal cavity was reduced, which influences the perception of nasal patency. It is suggested that wearers of high-efficiency masks that have minimal porosity and low air exchange for CO2 regulation should consider the amount of time they wear the mask.

Thermodynamic analysis and optimization of a high temperature triple absorption heat transformer.

First law of thermodynamics has been used to analyze and optimize inclusively the performance of a triple absorption heat transformer operating with LiBr/H2O as the working pair. A thermodynamic model was developed in EES (engineering equation solver) to estimate the performance of the system in terms of the most essential parameters. The assumed parameters are the temperature of the main components, weak and strong solutions, economizers' efficiencies, and bypass ratios. The whole cycle is optimized by EES software from the viewpoint of maximizing the COP via applying the direct search method. The optimization results showed that the COP of 0.2491 is reachable by the proposed cycle.